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Abstract 

There is evidence that the transiting planet HD 209458b has a large exosphere of neutral hydrogen, 
based on a 15% decrement in Lyman-a flux that was observed by Vidal-Madjar et al. during transits. Here 
we report upper limits on Ha absorption by the exosphere. The results are based on optical spectra of 
the parent star obtained with the Subaru High Dispersion Spectrograph. Comparison of the spectra taken 
inside and outside of transit reveals no exospheric Ha signal greater than 0.1% within a 5.1 A band (chosen 
to have the same AA/A as the 15% Lya absorption). The corresponding limit on the column density of 
n = 2 neutral hydrogen is N2 10 9 cm~ 2 . This limit constrains proposed models involving a hot (~10 4 K) 
and hydrodynamically escaping exosphere. 
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1. Introduction 

A milestone in extrasolar planet research was reached 
when Charbonneau et al. (2000) and Henry et al. (2000) 
observed the photometric signal of transits by the low- 
mass companion of HD 209458. The companion was orig- 
inally discovered by radial velocity measurements (Mazch 
et al. 2000), which specified only the orbital period (3.5 
days), orbital eccentricity (e < 0.03), and minimum mass 
of the companion (Ms'mi = 0.69 Mj up ). The transit light 
curves allowed the measurement of the companion's mass 
(by breaking the sini degeneracy) and radius (from the 
depth of the transit), providing an unambiguous case of a 
planetary-mass companion to a main sequence star, and 
a demonstration that the companion's density was that of 
a gas giant planet. 

Aside from this historic importance, the discovery made 
possible a number of unique and important follow-up 
studies, based on more subtle changes in the received 
starlight that should occur during transits. In this pa- 
per we are concerned with changes due to the passage of a 
small fraction of the starlight through the planetary atmo- 
sphere. Selective absorption by atmospheric constituents 
causes the transit depth to depend upon wavelength, an 
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effect that forms the basis of "transmission spectroscopy." 
Upper limits on various atmospheric absorption features 
have been given by Bundy and Marcy (2000), Moutou et 
al. (2001), Brown, Libbrecht, and Charbonneau (2002), 
and Moutou et al. (2003). The first successful detection 
using this technique was by Charbonneau et al. (2002), 
who observed a (0.023± 0.006)% increase in transit depth 
in the yellow light of the sodium resonance doublet. A 
strong sodium signal had been predicted by Seager and 
Sasselov (2000), and subsequently by Brown (2001) and 
Hubbard et al. (2001). 

More recently, Vidal-Madjar et al. (2003) reported a 
remarkably large transit depth of (15±4)% in the ultravi- 
olet light of the Lyman-a (Lya) transition of neutral hy- 
drogen. To produce such a large absorption, a spherical 
cloud of neutral hydrogen would need a radius of 4.3i?,j up 
or greater (depending on the optical depth), as compared 
to the planetary radius of 1.3i?j up inferred from broad- 
band observations, and the Roche lobe radius of 3.6i?j up . 
Such a signal could be produced by an extended enve- 
lope of hydrogen atoms that are evaporating from the 
planet's upper atmosphere. Vidal-Madjar et al. (2003) 
noted that the apparent blucshift of the absorbing atoms 
(up to —130 km s _1 ) provided additional evidence for at- 
mospheric escape, and speculated that this evaporation 
process accounts for the rarity of extrasolar planets with 
orbital periods smaller than 3 days. Additional data ob- 
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tained by Vidal-Madjar et al. (2004) were consistent with 
the previous Lya result, and provided evidence (at the 
2-3(7 level) for exospheric absorption by oxygen (O i) and 
carbon (C II). 

Further confirmation of the existence of this "exo- 
sphere," and characterization of its temperature, den- 
sity, and composition are clearly important goals. 
Unfortunately, theoretical predictions for the properties 
of the exosphere are not robust, making it difficult to de- 
sign strategies for observing the exosphere. One approach 
is to search for additional effects of hydrogen, which has 
already been implicated by the Lya measurement. If a 
significant fraction of the hydrogen exists in the first ex- 
cited state (n — 2), it will produce extra absorption in the 
Balmer lines of the stellar spectrum, and in particular Ha 
(6563 A). 

Thus we were motivated to search for Ha absorption 
from the exosphere of HD 209458b. A positive detection 
would confirm the existence of the hydrogen exosphere 
and constrain its density and temperature, with implica- 
tions for the long-term evolution of "hot Jupiter" extraso- 
lar planets. It would also represent a new and more prac- 
tical means of searching for exospheres of other extrasolar 
planets, given the necessity of observing Lya absorption 
above the Earth's atmosphere, and the confusing effects 
of Lya absorption by interstellar hydrogen and emission 
from the geocorona. Previously, Bundy and Marcy (2000) 
placed upper limits of 3-4% on H/3 and H7 absorption in 
a 0.3 A band, although they did not examine Ha due to 
contamination of that region of the spectrum by I2 lines. 
Moutou et al. (2001) also searched for exospheric absorp- 
tion in optical spectra. They did not remark on Ha in 
particular, but set upper limits of «1% for any features of 
width 0.2 A, and noted that the limits could be improved 
significantly with higher quality spectra taken in better 
atmospheric conditions. 

We have been obtaining optical echelle spectra of 
HD 209458 with the Subaru Telescope, over a wide range 
of planetary orbital phases, in order to search for reflected 
light from the planet. Data taken on the night of a plan- 
etary transit, described in section 2. are well-suited for 
the Ha search. The procedure by which we compared the 
spectra taken at different times is given in section 3, and 
significantly improved upper limits on exospheric Ha ab- 
sorption are derived from the data. In section 4 we discuss 
the physical interpretation and provide a brief summary. 

2. Observations and Data Reduction 

We observed HD 209458 on UT 2002 October 25 with 
the Subaru 8.2 m telescope of the National Astronomical 
Observatory of Japan, on Mauna Kca, Hawaii. The planet 
was predicted to transit its parent star in the middle 
of that night, according to the ephemeris of Brown et 
al. (2001). We used the High Dispersion Spectrograph 
(HDS; Noguchi et al. 2002) mounted at the //12.5 Optical 
Nasmyth focus. The entrance slit was 4" long and 
0'.'8 wide, oriented at a constant position angle. We 
used the standard Yb setup, in which the beam passes 



through an order-blocking KV370 filter, is collimated by 
the red-optimized mirror, dispersed by the echelle grating 
(31.6 grooves mm" 1 , blaze angle 71.°25), cross-dispersed 
by the red-optimized grating (250 grooves mm -1 , blaze 
angle 5.° 00), and detected on two CCDs, each having 
4100 x 2048 pixels of size 13.5 fim. The wavelength cover- 
age was 4100 A < A < 6800 A with resolution R « 90000. 

We obtained 29 spectra of HD 209458, each with an 
exposure time of 500 seconds, through air masses rang- 
ing from 1.0 to 1.9. Additionally, at both the start and 
end of this series of exposures, we obtained a spectrum 
with an I2 cell behind the slit (and a narrower slit width 
of 0'.'4) in order to determine accurate radial velocities 
for the star and verify the ephemeris. At the end of the 
night, we obtained 7 spectra of the B5 Vn star HD 42545 
in order to estimate the telluric spectrum: any narrow 
absorption features in the spectrum of this rapidly rotat- 
ing star (fsini = 245 km s ; Abt et al. 2002) are almost 
certainly telluric. The weather conditions were excellent, 
and the seeing was approximately 0"6. 

We used standard IRAF 1 procedures to process the 
frames and extract one-dimensional spectra. The frames 
were debiased, trimmed, and multiplied by the gain. A 
median flat field was created from 20 dome-lamp expo- 
sures and normalized in two dimensions (along the disper- 
sion, to take out the lamp spectrum, and perpendicular 
to the dispersion, to take out the aperture profile), and 
then used to correct for pixel-to-pixel sensitivity varia- 
tions. A two-dimensional smooth polynomial surface was 
fitted to the inter-order scattered light in each frame, and 
subtracted. The wavelength solution was determined from 
thorium-argon arc exposures taken at the start and end 
of each night. A total of 54 orders were extracted, after 
fitting the trace of each order of each spectrum separately, 
and summing the counts across the dispersion with opti- 
mal weighting. In the order spanning Ha, the pixel scale 
was « 0.02 A pixel -1 and the full-width at half-maximum 
of arc lines was approximately 4 pixels. The resulting one- 
dimensional spectra had a typical signal-to-noise ratio of 
« 300 in each pixel. 

Exposures taken with the I2 cell were processed sepa- 
rately, using the standard procedures described by Butler 
et al. (1996) and Sato et al. (2002). We determined the ra- 
dial velocity based on the blue CCD only. Figure 1 shows 
our measurements, which are uncertain by sa 6 m s , 
as estimated from the variance in results from 100 sepa- 
rate spectral segments in each exposure. The two points 
from October 25, obtained at the start and end of the 
HD 209458 observations on that night, bracket zero ra- 
dial velocity (orbital phase 0.5), confirming our expecta- 
tion that the transit took place between those times. 

1 The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the U.S. National Optical Astronomy Observatories, 
which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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3. Transmission Spectroscopy 

During the transit, the planet blocks a portion of the 
stellar surface, causing the flux received from the star to 
decrease. Our goal was not to measure the transit light 
curve, which would have required observations of a com- 
parison star to correct for flux variations due to Earth's 
atmosphere. Rather, our goal was to measure differences 
between the light curve observed in Ha and at other wave- 
lengths, effectively using part of the parent star's spec- 
trum as the comparison signal. For exosphcric effects to 
be detectable, they must not only be large enough to over- 
come the Poisson noise in the starlight, but they must also 
be distinguished from spectral variations caused by insta- 
bilities in the Earth's atmosphere and the spectrograph. 
In this section we describe the method we used to attempt 
to make this distinction. 

3.1. Corrections for Instrumental Variations 

Figure 2 compares two representative spectra that were 
taken 2.5 hours apart. The top three panels show three 
consecutive orders bracketing Ha, the strong absorption 
feature at 6563 A. We refer to these spectra as Si(X n j) 
and ^(An.j), where X n ,j is the jth wavelength bin within 
the nth spectral order. Note that we have not normalized 
the spectra; the continua are peaked near the middle of 
each order due to the blaze function of the spectrograph. 

The green points in the bottom three panels of figure 2 
show the ratios R(X n ,j) = Si(X n ,j)/S2(X n ,j)- The mean 
of the ratio spectrum is not unity, which is expected, be- 
cause of variations in total flux received in each exposure 
caused by the Earth's atmosphere. More surprising are 
the wavelength variations in the ratio spectra: there are 
±5% variations varying smoothly over sa 5 A scales. We 
believe these variations are not an artifact of the extrac- 
tion procedure, because the same variations were obtained 
from spectra that were extracted using different aperture 
sizes, and from spectra that were extracted using custom- 
built procedures in the Interactive Data Language (IDL) 
rather than IRAF. 

The variations are surely instrumental, because the pat- 
tern is very similar in each order: the ratio spectrum is a 
function of j rather than A, and depends on n to a much 
lesser degree. The ratio spectra also exhibit a smooth time 
dependence. In particular, the wavelength variations are 
greatest when one of the spectra involved in the ratio was 
taken at the beginning of the night. This was when the 
target was near zenith and the instrument rotator was 
moving most quickly. The ratio spectra are flatter when 
the two spectra were taken with nearly the same rotation 
angle. For these reasons, we hypothesize that these appar- 
ent variations in the blaze function are caused by flexure 
of the spectrograph. Similar variations have been seen 
by other observers using Subaru/HDS (e.g., Nagao et al. 
2003) and also Keck/HIRES (e.g., Suzuki et al. 2003). 

We corrected for this effect empirically by using the pat- 
tern observed in the adjacent orders. To correct the nth 
order, we computed i?(A„+i.j) and R(X n -ij), and then 
smoothed each of these functions in j with a boxcar size 



of 100 pixels (2 A), giving R{n + l,j) and R{n — l,j)- We 
then used linear regression to determine the two numbers 
c„ + i and c„_i that minimized the sum of squared residu- 
als between S\ and 

R(n+l,j)+c n - l R(n~l,j)] 5 2 (A nii ).(l) 

For a given pair of spectra we typically found c„+i ~ 
Cn-i ~ 0.5. The resulting "blaze corrected" ratio spec- 
trum, Si/iS^bj is near unity with a standard deviation 
within 10% of that expected from the Poisson statistics 
(see black dots in figure 2). 

Next, we corrected for small variations in the wave- 
length scale by determining the two numbers cq and AA 
that provided the best match between Si and 



5 , 2m(A„ J ) — C S2b(X n j) + A\ 



dS 2h (X) 
dX 



(2) 



where the derivative was approximated numerically by 
3-point Lagrangian interpolation. The wavelength shifts 
were always smaller than 0.005 A (0.25 pixel). We refer 
to S^m as the "matched" spectrum; it has been matched 
in blaze function and wavelength scale to Si. 

3. 2. Calculation of Difference Spectra 

Having established this matching procedure to correct 
for instrumental variations between any pair of spectra, 
we created a template spectrum with a high signal-to- 
noise ratio in the following manner. A spectrum from the 
middle of the night was chosen as a reference spectrum. 
All the spectra taken through an air mass less than 1.4 
were matched to the reference spectrum. (Spectra taken 
through larger air masses were significantly contaminated 
by time- variable telluric lines, as will be seen shortly.) The 
preliminary template spectrum was defined as the median 
of all the matched spectra. Then, the entire procedure 
was repeated, using the preliminary template spectrum 
as the reference spectrum, resulting in the final template 
spectrum T(X). Finally, a time series of difference spec- 
tra was created, by subtracting each spectrum from the 
suitably matched template. 

Some examples of the difference spectra are shown in 
figure 3, for a 55 A region in the vicinity of Ha. The 
difference spectra are plotted in units of Poisson deviates, 

S(X)-T m (X) 



with the space between minor tick marks indicating one 
Poisson deviation. Also plotted are the template spec- 
trum (bottom) and the estimated telluric spectrum (top). 
To produce the telluric spectrum, we created a template 
spectrum for HD 42545 in the same way as we did for 
HD 209458. Then we removed the rotationally broad- 
ened Ha absorption line (which had a full-width at half- 
maximum of ~11 A) by fitting a smooth polynomial to the 
spectrum and dividing by this polynomial. We note that 
the telluric spectrum is shown for visual reference only; we 
did not use it in any of the calculations described below. 

The difference spectra are nearly consistent with 
Poisson noise: they have standard deviations near unity, 
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as expected. However, a close inspection of figure 3 shows 
the residuals are not truly random. For spectra taken 
through large air masses (late times), telluric lines are vis- 
ible. In addition, there are excursions of order 1— 2a that 
are evident as "ripples" on scales of 2-5 A (100-250 pix- 
els). For example, the spectrum at t = 0.24 shows a saw- 
tooth pattern in the residuals between 6535 and 6545 A. 
We do not know the cause of these ripples, but we are con- 
fident that they are not due to the exosphere of the planet, 
for two reasons: they are not specific to any particular ab- 
sorption line in the stellar spectrum; and they also appear 
in spectra taken on October 27, when the planet was not 
transiting. 

3.3. Calculation of Difference Light Curves 

There are no obviously significant features in the dif- 
ference spectra at the position of Ha, implying that once 
the spectra have been adjusted to correct for total flux and 
instrumental effects, there is no significant difference be- 
tween the Ha light curve and those of neighboring bands. 
To plot these light curves, we computed the fractional dif- 
ference in flux between the spectrum and the matched 
template, 

f 2 d\ [S(X,t)-T m (X)} 

m = ^— ^ (4) 

/ d\ T m (X) 

for the 3 different wavelength bands indicated in figure 3: 
a band centered on Ha (green), and two comparison bands 
(red and blue) in the adjacent continua. The width of the 
Ha band (5.13 A) was chosen to match the AA/A of the 
band over which Vidal-Madjar et al. (2003) reported Lya 
absorption. The central wavelength of the Ha band was 
shifted to compensate for the time-variable Doppler shift 
of the planetary orbit. This shift is evident as a tilt in 
the green lines of Figure 3. The central wavelengths of 
the comparison bands were chosen to avoid telluric lines 
and strong stellar absorption lines, as much as possible. 
The widths of the red (3.87 A) and blue (3.61 A) com- 
parison bands were chosen such that all three bands have 
the same total flux in the template spectrum. Figure 4 
shows the light curve S(t) for each of the three bands. 
Apart from a few outliers, the variations in the Ha band 
are similar in size and character to the variations in the 
comparison bands. The red and blue light curves exhibit 
small (and opposite) gradients which are probably effects 
of increasing air mass. 

To place a quantitative limit on excess Ha absorption, 
we calculated <5i n , the mean value of 6(t) for spectra taken 
during transit (between second and third contacts). As 
an estimate of the uncertainty in 5- m , we computed the 
standard deviation of S(t) for the 12 in-transit spectra 
and divided by -\/l2- Likewise, we computed <5 ut and its 
uncertainty from the spectra taken before first contact or 
after fourth contact. 2 Then we computed 6n a = din — ^out, 
which would be negative if there were extra Ha absorption 

2 To compute the transit times, we used an ephemeris defined by 



during transit. For the 5.13 A band, the result is Sn a = 
-0.00024 ±0.00029. We investigated the dependence of 
this limit upon bandwidth, by repeating this procedure 
for bands varying in width from 0.5 A to 10 A, centered 
on Ha. The results are shown in the upper right panel of 
figure 4. Finally, we tried blueshifting the 5.13 A Ha band 
by —130 km s _1 , to match the possible blueshift in Lya 
observed by Vidal-Madjar et al. (2003). The outcome was 
also a null result, 5 Ha = 0.00015 ± 0.00021. 

Although our spectrum-matching procedure was in- 
tended to correct only for instrumental and telluric effects, 
it also acts to dilute any real signal. To estimate the size 
of this effect we injected an artificial signal into the data. 
We multiplied the extracted spectra by (1 — g\), where g\ 
is a Gaussian function centered on Ha with a full-width 
at half-maximum of 3 A (a fairly arbitrary choice, since 
the true width of the absorbing feature is unknown), rep- 
resenting a total flux decrement of 0.1% in the 5.13 A 
band over which Sn a was measured. When the preced- 
ing analysis was performed on the adulterated data, the 
result was 5 Ha = -0.00081 ± 0.00028. This is greater by 
0.00057 than the measurement with no artificial signal, 
verifying that the artificial signal (although diluted) could 
be recovered in spite of the noise. We conclude than an 
extra 0.1% absorption during transit would have been de- 
tected. Figure 5 shows the resulting light curve 5(t) and 
the "measured" absorption as a function of the width of 
the Ha band. 

4. Discussion and summary 

Fundamentally, our null result corresponds to a max- 
imum number of neutral hydrogen atoms in the n = 2 
state that are present in the exosphere of HD 209458b. 
Likewise, the detection of Lya absorption by Vidal- 
Madjar et al. (2003) implies a minimum number of neu- 
tral hydrogen atoms in the n = 1 state. Unfortunately 
it is difficult to use these results to constrain interesting 
quantities such as the exosphcric temperature and den- 
sity, because the physical conditions in the exosphere are 
likely to be quite complex. The column density of neutral 
hydrogen may vary widely across the surface of the pri- 
mary star due to evaporation, tidal forces, and radiation 
pressure (see, e.g., Lecavelier des Etangs et al. 2004). Gas 
may be present at a variety of temperatures and densi- 
ties, and may be far from local thermodynamic equilib- 
rium (LTE). [Even in the denser atmosphere, non-LTE 
effects may be required to explain the unexpectedly weak 
detection of sodium by Charbonneau et al. (2002), as pro- 
posed by Barman et al. (2002) and Fortncy et al. (2003).] 
Obviously the limited empirical information available is 
not sufficient by itself to determine a realistic physical ex- 
osphere models. 

With these caveats, we can use a simple order-of- 
magnitude argument to provide an upper limit on the 
column density of n = 2 neutral hydrogen that exosphere 

J.D. (mid-transit) = 2451659.93675 (Brown et al. 2001) and the 
most precise published period, P = 3.524739 d (Robichon and 
Arenou 2000). 
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models will need to obey. In order to translate our upper 
limit on additional Ha absorption into a corresponding 
upper limit on the Ha equivalent width, we combined the 
in-transit spectra to form 5; n (A), and we combined the 
out-of-transit spectra to form S out (A). Then the equiva- 
lent width W\ of any transit absorption feature is 



dX 



Sout (A) — S^n (A) 
Sout(A) 



(5) 



Performing the integral over the 5 A band described in 
the previous section, we find Wn a < 1.7 mA. If we further 
assume that the fraction of the stellar surface covered by 
the exosphere is AA/Am 0.15 (as estimated from the Lya 
result), then the null result implies that Ha absorption is 
weak, and the corresponding limit on the column density 
of n = 2 atoms (N2) is determined via 



AA\ ire 2 



/23Ah q -^2, 



(6) 



where /23 ~ 0.641 is the oscillator strength of the n = 2 — > 3 
transition. The resulting limit is N2 < 1.0 x 10 9 cm~ 2 . 

In principle, Ni could be estimated from the Lya decre- 
ment observed during transits, and the ratio N^/Ni could 
be related to the H I excitation temperature. In reality, 
Ni is highly uncertain because the absorption is probably 
saturated, and the line profile is unobscrvablc due to in- 
terstellar absorption and geocoronal emission. However, 
any exosphere model that specifics ATi is constrained by 
our results to have a maximum excitation temperature 
T cx , according to the relation 



A2 

JVi 



92 

I 

91 



■ cxp 



E2 — E\ 



4 exp 



10.2 eV 



(7) 



This constraint appears to be important in the context of 
recently proposed models in which the exosphere is very 
hot and hydrodynamically escaping, as argued below. 

The temperature of the lower atmosphere of 
HD 209458b is likely to be near the radiative effec- 
tive temperature, T e g « 10 3 K, at which the n = 2 
population predicted from equation (7) is utterly negli- 
gible. But, as argued by Moutou et al. (2001), Lammer 
et al. (2003), and Lecavelier des Etangs et al. (2004), the 
exosphere is likely to be significantly hotter than T e g, 
just as the exosphere of Jupiter is considerably hotter 
than its lower atmosphere (^10 3 K vs. 150 K; see, e.g., 
Atreya 1986). Lammer et al. (2003) found that close-in 
gas giant planets could have exospheric temperatures up 
to 10 4 K due to intense X-ray and extreme-ultraviolet 
(EUV) irradiation. Depending on the density of neutral 
hydrogen and other factors, such a hot exosphere could 
produce detectable Ha absorption in violation of our 
constraint. 

For example, Lecavelier des Etangs et al. (2004) pre- 
sented 3 models for the thermosphere and exosphere of 
HD 209458b. In each case they specified the density pro- 
file of neutral hydrogen and computed the size and shape 
of the exosphere as a function of kinetic temperature Tk, 
taking tidal distortion into account. For the specific case 
of Model A, the density is = 2 x 10 9 cm -3 at the 



base of the thermosphere, and decreases with radius ac- 
cording to the barometric law. We computed the aver- 
age column density Ni within the annulus reaching from 
the thcrmobasc to the Roche radius, and then used equa- 
tion (7) to translate our upper limit on N2 into an upper 
limit on the excitation temperature, finding T ox < 8000 K. 
By comparison, Lecavelier des Etangs et al. (2004) found 
that Tk > 8000 K in order to produce an escape rate 
of 10 10 g s _1 [the lower limit inferred by Vidal-Madjar 
et al. (2003)], and T k > 11100 K through a simulation 
of EUV heating. Of course, one would not necessarily 
expect T cx w Tk in the tenuous thermosphere and exo- 
sphere, because of departures from LTE; in general, T cx 
is neither an upper bound nor a lower bound on Tk- The 
n = 2 state may be depopulated when collisions are infre- 
quent, and likewise, whatever process accelerates the gas 
to ~100 km s^ 1 may overpopulate the n = 2 state via re- 
combination cascades and Lya resonant trapping. A full 
non-LTE exosphere model will need to accommodate the 
resulting upper limit on the excitation temperature. 

In summary, we have used high resolution and high 
signal-to-noise ratio Subaru-HDS optical spectra to search 
for excess Balmer Ha absorption during a transit of the 
extrasolar planet HD 209458b. No excess absorption was 
detected. Our upper limit is two orders of magnitude be- 
low the Lya absorption reported by Vidal-Madjar et al. 
(2003, 2004). It may be difficult to improve upon this 
limit with ground-based instruments, given the difficulty 
of correcting for telluric and instrumental spectral varia- 
tions. However, there should be no obstacle to improving 
upon our limit with space-based spectroscopy, and we are 
aware that such an effort is underway using data from 
the Hubble Space Telescope (Charbonneau, D., private 
communication). The current result, and any future re- 
finements, will be useful in the further development of 
models for the escape of the hot exospheres of close-in gas 
giant planets. 

We are very grateful to Yutaka Abe, Dave 
Charbonneau, Bruce Draine, Bob Kurucz, Avi Loeb, 
Paul Martini, Dimitar Sasselov, and Motohide Tamura, 
for helpful consultations. This work is based on data 
from the Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan. We wish to 
recognize and acknowledge the very significant cultural 
role and reverence that the summit of Mauna Kea has 
always had within the indigenous Hawaiian community. 
We are most fortunate to have the opportunity to conduct 
observations from this mountain. 
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Fig. 1. Subaru radial velocity measurements of HD 209458 on 2002 October 25 and 27. The full radial velocity curve is also 
plotted, using the orbital parameters determined by Naef et al. (2004). This confirms that the planetary transit took place be- 
tween the times of our two measurements on October 25, which bracketed the time series of spectra presented in this paper. 




Fig. 2. The top panels show sample spectra of HD 209458 taken 2.5 hours apart. Each panel shows one order. 
The bottom panels show the ratio of the spectra, both before (green) and after (black) the blaze function correction. 
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Fig. 3. Differences between individual spectra of HD 209458 and the matched template, expressed in units of Poisson deviates. The 
spacing between minor tick marks represents one Poisson deviation. For clarity, only 14 of the 29 spectra are shown. The template 
spectrum is depicted at the bottom of the series, and the telluric spectrum is depicted at the top. The right-hand axis shows the time 
each spectrum was taken, in the same units as figure 4. Colored lines mark the edges of the three band passes described in the text. 
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Fig. 4. Comparison of the Ha light curve and light curves in adjacent bands. (Top left) The fractional difference in flux between each 
spectrum 5 and the corresponding matched template T m , within the 5.1 A Ha band shown in figure 2. Solid lines indicate the times of 
first and fourth contact. Dashed lines indicate the times of second and third contact. (Bottom left) Same, but for the red comparison 
band. (Bottom right) Same, but for the blue comparison band. (Upper right) The difference between in-transit and out-of-transit 
flux deviations from the template, as a function of the width of the band pass. The single filled point is for the 5.1 A band. 
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Fig. 5. Same as the top two panels of figure 4, except that an artificial absorption signal of strength 0.1% has been added to the 
in-transit spectra. From this test we concluded that a 0.1% excess in Ha absorption during transits would have been detected. 



